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Abstract. Higher-rank Minkowski valuations are efficient means for describing the 
geometry and connectivity of spatial patterns. We show how to extend the framework of 
the scalar Minkowski valuations to vector- and tensor- valued measures. The versatility 
of these measures is demonstrated by using simple toy models as well as real data. 



1.1 Introduction 



The spatial patterns originating from the polymorphic aggregation of matter in 
Nature occur on vastly different length scales and with unlimited variety. For 
instance, the pattern shown in Figure |l.l| could represent a biological structure 
obtained with X-ray tomography, or a colloidal cluster in a milky emulsion. In 



fact, Figure 1.1 displays an isodensity contour of cosmic matter as traced by the 
galaxies surrounding the Milky way, which lies in the center of the plot. 
The visual impression from such pictures may already convey valuable insight 
into the dynamic evolution or - in the context of biology - into the functionality 
of such structures. However, for unbiased inferences as well as for the compar- 
ison with model simulations, one needs objective and quantitative measures to 
characterize the geometry and topology of typical structural motifs. 
In this article we deal with a family of functions which provide a unique descrip- 
tion of the morphology displayed by spatial patterns. We call these functions the 
Minkowski valuations; they are extensions of the well-known Minkowski func- 
tionals often mentioned within this volume (see the contributions by C. Arns et 
al. and M. Losche and P. Kriiger in this volume as well as |2f| in the previous 
volume LNP 554). Therefore they share some of the simplicity and beauty which 
the Minkowski functionals own. In particular, they are characterized by simple 
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Fig. 1.1. The matter distribution around our galaxy as traced by the PSCz survey. 
The galaxy numbe r density is smoothed; one isodensity contour is shown (from the 
PSCz homepage at http : / / star- www, dur. ac , uk / cosmology / theory /pscz , htmj ) . 



axioms which set intuitively reasonable standards for a morphological descrip- 
tion. In the same way as the Minkowski functionals generalize the notion of the 
volume, the Minkowski valuations extend the notion of the mass vector (or cen- 
ter of mass) and of the inertia tensor. Like their scalar relatives the majority of 
the Minkowski valuations are curvature measures integrated over the surface of 
a pattern. 

We start with a short motivation in Section |1.2j by showing that there are good 
reasons to go beyond the scalar Minkowski functionals. After having outlined 
ways to generalize the Minkowski functionals (Section 1.3), we review some of 



the extensions' mathematical properties in Section 1.4. To illustrate how useful 



the Minkowski valuations are, we discuss concrete applications covering sub- 
molecular up to cosmological scales (Sections 1.5 and 



1. 



Before, two cautionary remarks are appropriate: first, in order to distinguish 
conceptually between the scalar Minkowski functionals and their higher-rank ex- 
tensions, we shall address the whole entirety of Minkowski functionals and their 
extensions via the term "Minkowski valuations" (MVs), whereas "Minkowski 
functionals" (MFs) is reserved to the scalar Minkowski functionals as before. 
Secondly, there is a vivid theoretical development of higher-rank Minkowski val- 
uations these days Therefore, this overview can not provide a complete 
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survey on Minkowski valuations; on the contrary, substantial progress is to be 
expected both from the mathematical side as well as from applications. For 
instance, morphological analysis based on Minkowski valuations may be use- 
ful to quantify molecular orientations (see the contributions by F. Schmid and 
N. H. Phuong, U. S. Schwarz and G. Gomppcr, D. Jculin, C. Arns et al., and 
H.-J. Vogel in this volume). Our vector-valued descriptors may also be used to 
refine the morphological analysis of Langmuir monolayer phases (presented by 
M. Losche and P. Kriiger in this volume) which has so far been based merely 
on scalar MFs. Another application may cover the geometric characterization of 
liquid foams described by F. Graner in this volume, where stress tensors need 
to be related to a complex spatial structure. 



1.2 Going beyond scalar Minkowski functionals — a 
physical motivation 



The following problem is typical for challenges occurring in physics: the forma- 
tion of galaxy clusters depends in a complicated way on the background cosmol- 
ogy usually described in terms of a cosmological model. Such a model is defined 
through the values of the cosmological parameters and the power spectrum of 
the initial density perturbations |30|H . In order to investigate the impact of the 
background cosmology on galaxy clusters one carries out gravitational TV-body 
simulations; typical results comprise cluster images such as those shown in Fig- 
To compare the models and to check which model fits best present-day 



1.2 



ure 

cluster observations, one has to use quantitative and discriminative morpholog- 
ical measures. 

The Minkowski functionals certainly place practical tools at the physicist's dis- 
posal. Defined on the convex ring /C^, i.e. for finite unions of compact, convex sets 
in ci-dimensional Euclidean space, they are the only motion-invariant, convex- 
continuous, and additive descriptors. Despite this very general characterization 
there are only a finite number of Minkowski functionals, namely (d + 1) in d 
dimensions (characterization theorem, [ fL6U2l|| ) . They have simple geometrical 
meanings derived from their representation as integrals. These integrals either 
extend over the pattern itself to yield the volume, or cover its surface which 
is weighted with combinations of the local curvatures in order to give - in two 
dimensions, e.g. - the perimeter and the integrated curvature (which equals the 
Euler characteristic). The Minkowski functionals also arise quite naturally as 
expansion coefficients in the Steiner formula |42f specifying the volume of the 
parallel body of a convex body. Indeed, the Minkowski functionals have been 
extensively applied to physical structures at vastly different length scales, see 
[ pH| for an overview and pp[ | for astrophysical applications. On the mathemati- 
cal side, numerous useful results such as the principal kinematic formulae or the 
Crofton formulae (see (l(|) have been proven with integral geometric methods. 
However, the morphometry with Minkowski functionals is by no means compre- 
hensive. Combining two of their defining principles, additivity and homogeneity, 
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Fig. 1.2. The projected matter distribution of two simulated galaxy clusters (from the 
GIF-project). Their different morphologies are due to different values of the cosmolog- 
ical parameters and to different initial power spectra imposed (left panel: r Cold Dark 
Matter model; right panel open Cold Dark Matter model, see |g{ for more details). 
Apart from that, however, the clusters stem from comparable initial conditions. In or- 
der to distinguish between both models we need a thorough morphological description. 



one can easily see, that, e.g., moving a connected part of a pattern results in sig- 
nificant visible changes which can not be detected by the Minkowski functionals 
at all. This is illustrated in the first panel of Figure 1.3 representing schematic 
clusters. It is surely of interest, however, in which way the subclumps of a system 
are arranged, and how they are aligned with respect to each other (see the second 
panel of Figure 1.3). Therefore, we have to move beyond the scalar Minkowski 
functionals. It turns out that one can generalize the Minkowski functionals with- 
out leaving the mathematical framework of integral geometry. 



1.3 The hierarchy of Minkowski valuations — extending 
the framework 

There are several ways to extend the framework of the Minkowski functionals. 
For instance, one may start with Steiner's formula. 

The Steiner formula. It is a natural choice to employ the volume V(K) = J„ dV 
of a body K for its very coarse description. If the parallel body K t is considered, 
which contains all points closer to K than e > 0, and if K is convex, then one 
finds from Steiner's formula an expansion for V(K £ ) in powers of e, 
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Fig. 1.3. Neither the schematic clusters displayed in panel l.a. and l.b. nor both 
patterns in panel 2. a. and 2.b. can be distinguished from each other using only the 
scalar Minkowski functionals. Additivity and motion-invariance prevent the MFs from 
discriminating between the different positions (upper panels) or different orientations 
(lower panels) of the subcomponents, respectively. 



where the Minkowski functionals W v arise as expansion coefficients. The latter 
functionals are useful morphometric concepts and can be extended to non-convex 
patterns. In this sense, the volume quite generally leads to a family of related 
measures. 



Since we saw that in some applications (recall the upper panels of Figure 1.3) 
the position of a body may be worth considering, we specify the location of K 
with its mass vector 



V 1 (K) = / dVx . (1.2) 
Jk 

Alternatively, we could introduce its center of mass po = V 1 /V . 



Moreover, as the lower panels of Figure 1.3 indicate, the body's orientation may 



be relevant. In this case we have to move on to even higher moments; for instance, 
the tensor 

V 2 {K)= [ dVxx (1.3) 

(where xx denotes the symmetric tensor product^) contains information on main 
directions of the body K. Note, that both the mass vector of a body and its 

1 As H show, there are no non-trivial antisymmetric Minkowski valuations of rank 
r > 1 at all, i.e. no motion-covariant, addit ive, and conditionally continuous tensor 



functions with rank r > 1, see below Section 1.3 for the explanation of these concepts. 
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inertia tensor (the traceless part of Equation [T^: V 2 (K ) — Tr (V 2 (K)} Ed, where 
Ed is the d-dimensional second-rank unit tensor) are of major interest for the 
physicist: the mass center of a rigid body moves in space as if the whole force 
acting on the body was concentrated onto a point particle at its center of mass. 
The inertia tensor determines how a body is rotating, given a fixed angular 
momentum. - Higher-rank moments such as symmetric r-rank tensors of the 
form V r (K) = J dVx r with r > 2, may be considered, too. 
Now it is again possible to write down a Steiner-type formula for the transition 



to the parallel body of a convex set K |40||43(] : the scalar expansion coefficients 



have to be replaced by vector- or tensor-valued coefficients. The following simple 
argument shows that not only the character, but also the number of coefficients 
has to change: for a unit ball centered at the origin, -Br(O) = {x||x| < i?}, the 
construction of the parallel body results in B( £+R ) (0) . On the other hand, as 
a simple calculation shows, V r (S( e+ m(x)) cx (e + R) d+r . Thus, powers up to 
(d + r) arise in the generalization of Steiner's formula, which we write formally 
as 

d+r 

V r (K e ) = J2 e " M ^ ■ (1-4) 

As for the scalar Minkowski functionals, one can obtain an integral representation 
of the Stcincr coefficients by considering a smooth and convex body K. From 




Fig. 1.4. The fact that each volume element in K is weighted with the position vector 
in order to give V r (K) for r > 0, makes the local normal n relevant for the parallel 
body K t . 



the integral representation of the scalar MFs and from the form of the moments 
V r one might guess that the coefficients M rv are built by functions of the form: 

dS^^Su (/Sl, .., Kd-l) (1.5) 

dK 

with the elementary symmetric functions s v for v = 0, .., {d — 1) of the local 
curvatures k\, .., Kd-i (for a definition, see Q). But this will not suffice, since 



we need (d + r) coefficients. A closer look shows (see Figure 1.4), that indeed 
the local normal vector comes into play; for, if x C dK, then x + en lies on 
the surface of K e ; here n denotes the local normal vector to dK at the point 
x. The weighting with the local position in the definition of V r therefore leads 
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to contributions from local normal vectors in the Steiner formula. In [Q it is 
shown that in general the expansion coefficients M TV can be written as linear 
combinations of the quantities: 

Wr = ~Lf dS d - 1 s^ 1 (K 1 ,.., Kd _ 1 )xPn" (1.6) 

V U JdK 

for v = l,..,d, and p,q € No- Now the W^ q look more basic than the coefficients 
arising in the Steiner formula for the volume tensor. It seems therefore, as if, 
in contrast to the case of the scalars, the Steiner coefficients are no longer the 
basic quantities. In order to generalize the Minkowski functionals we therefore 
have to move to an alternative route. 



Spatial moments of the Minkowski functionals. Let us take the moments V r 
together with the W£' s (u = 0, .., (d — 1), r, s = 0, 1, ..) as basic quantities. In 
order to arrive at a coherent notation, we set Wq° = V r and Wq ' j = for j > 0. 
Higher-rank Minkowski functionals therefore are higher moments of the scalar 
Minkowski functionals where the local surface is weighted with both position 
and surface normal. If the body K is shifted with a translation vector t, then 
the WJ'; S transform like 




as follows immediately from their definition. For rotations R, we have 

Wl' s (RK) = RW^ S (K) , (1.8) 

where we use the fact that the <i-dimensional orthogonal group 0(d) has a natural 
representation in each tensor space. 

Recalling that the Minkowski functionals can be described axiomatically, one 
would be inclined at this point to look for a characterization theorem for vectors 
and tensors. It turns out that one has to include even more valuations in order 
to achieve that purpose. 



Axiomatic characterization. A rigorous way to extend the Minkowski functionals 
is therefore to relax one of their defining axioms. We keep the additivity 

<P(A U B) = <P(A) + <P(B) — <P(A n B) (1.9) 

for a tensor valuation we also maintain the conditional continuity, i.e. continu- 
ity for the subclass of convex bodies^; but, in the spirit of our previous findings, 
we replace the motion invariance by motion covariance defined as follows: for 

2 On the space of convex bodies one uses the Hausdorff metric, see p^ ], p. 48 for a 
definition. 
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translations we require the existence of a number n 6 No and of functions <Pi 
(i = 0, .., n) such that for all K 

n 

${K + t) =^T$ i (K)t i . (1.10) 

i=0 

Note, that products such as <Pi(K)V always denote symmetric tensor products. 
For rotations we use the natural representation of 0(d) in the space of tensors 
as before, and require for R G 0(d): 

$(RK)=m(K) . (1.11) 

In a seminal work, Alesker (J] proved the following characterization theorem: ev- 
ery motion-covariant, convex-continuous, and additive mapping <P can be written 
as a linear combination of the functionals 

E k d W r „> s (1.12) 

with i, j, k G No and v = 0, .., d. Ed denotes the d-dimensional unit matrix^]. 
Equation ( jL12j ) gives all of the Minkowski valuations and finally answers the 
question how to generalize the MFs in a suitable manner. It contains more 



valuations than those given in Equation (1.6). However, some caution has to 
be exercised: Alesker's theorem Q only bounds the dimension of the tensor 
functional spaces from above. It is not obvious that all of the above functionals 
are linear independent. Indeed, pif proved some linear relationships holding 
among the tensors. In order to give one example, Gaufi' theorem yields 

W 1 ' 1 = I dS d - 1 xn = W°' t) (K)E d . (1.13) 

JdK 

One therefore has to explore the different tensor ranks further. 



1.4 Exploring higher-rank Minkowski valuations 



1.4.1 Some mathematical results 



Minkowski vectors. Minkowski vectors, i.e. Minkowski valuations of rank one, 
have been known for some years p^ , ^p| , p]p^ ] . There is a strict one-to-one cor- 
respondence between the vectors and the scalars: for each Minkowski functional 
W v oc J gK s„_i (k\, .., Kd-i) there is a related vector j gK s„_i (ki, .., Kd-i) x, 
such that (d+1) vectors exist altogether. The reason why no normals enter at this 
level is basically because for compact convex bodies J gK s v -\ (k\, .., k^-i) n = 

3 For higher than two-dimensional spaces, it does not make a difference whether we 
require 0(d)- or S'0(d)-covariance. In the two-dimensional case, however, it does ^J. 
For the one- dimensional case see also H. 
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for v = 0, .., (d — 1) (for an elementary proof s ec fll7| |). The integral geometry of 
vectors was extensively scrutinized, see, e.g., ]18[]. For physical applications of 
vectors (see p| Jlf|) the curvature centroids 



Pi = W^ /W^° (i/ = 0,..,d) (1.14) 

are a useful concept. Their meaning is straight forward: they locate morpho- 
logical information. We shall see below that, if the curvature centroids do not 
coincide within one point, strong evidence for asymmetry is given. 



Second-rank tensors. The situation becomes more complicated for higher-rank 
tensors. It can be shown, that at the second rank, the following (3c? + 1) tensors 
provide a base [|] : 

d^xx , 

/v 

dS^-^xx , f dS^SuXn , f dS^s^nn (1.15) 

dK JdK JdK 

for v = 0, .., (<£— 1). Not all of them are interesting from a physical point of view, 
since some of them are related to Minkowski functionals times the unit tensor 
and thus do not contain any new information. Principal kinematic formulae and 
Crofton formulae are known for second-rank tensors p3j,pp3] . 
For a geometrical interpretation, we have to bear in mind, that the values of the 
tensors depend on the origin adopted. For W%'° (y — 0, .., d), however, the cor- 
responding curvature centroid p v is a natural reference point; the other tensors 



listed in Equation (1.15) are translation-invariant. Obviously, with such a choice 
the tensors are sensitive to morphological anisotropies, in the case of symmetry 
they align along the symmetry axes, otherwise they distinguish morphologically 
interesting directions. It is therefore useful to consider their eigenvalues and 
eigendirections . 



Beyond second rank. Although higher-rank tensors may be interesting as well, 
we shall not consider them further at this point. 



1.4.2 Examples 

In order to give a first example of how the Minkowski valuations discri min ate 
spatial structures, we consider a few simple patterns as shown in Figure |l.5| . In 
the case of a circle (panel 1) the curvature centroids coincide within the circle's 
center; the tensors - here represented as an ellipse featuring their eigenvalues 
and the corresponding directions - are isotropic. It is useful to divide the ten- 
sors W^' s by the corresponding Minkowski functionals W®'°. The trace of the 
normalized tensor (i.e. the sum of the eigenvalues) quantifies to which extent 
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the morphological information is concentrated around the center, how large the 
morphological fluctuations around the center are. 

The transition to an ellipse of the same volume as the circle (panel 2) can be 
traced either with the tensors becoming anisotropic or with the ratio of the 
squared surface to the volume (isoperimetric ratio, see J3qj ). Furthermore, as 
the morphological concentration is decreasing, the traces of most of the normal- 
ized tensors increase, especially of W 2 : °: the more elongated the ellipse becomes, 
the higher is its curvature at its very fringes along the horizontal direction. The 
normalized trace of W-y , on the contrary, becomes smaller, since the normal 
and the local position vector cease to be parallel to each other. The curvature 
centroids still coincide with each other because of the point symmetry. - If we 
reduce the symmetry further to a merely axially symmetric configuration (panel 
3), the curvature centroids fan out along the symmetry axis. In the case of no 
symmetry at all (panel 4) they may constitute a triangle. 

Finally, suppose one is given three convex bodies (or d + 1 ones in d dimensions) 
with known Minkowski functionals in order to combine them into a pattern 
such as indicated in panels 5 and 6. Given the scalar measures of the whole 
pattern, especially its Euler characteristic, one can constrain the number of in- 
tersections or subclumps in a pattern, which does not make any difference in 
our examples. For disentangling the positions of the grains, however, one has to 
know the curvature centroids. The centroids weight the subclusters according 
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to their Minkowski functionals. Therefore, whereas from the volume's point of 
view, only the larger components contribute significantly, all subclumps have 
the same weight (Euler characteristic = 1) for the curvature centroid P2. If the 
subclumps do not overlap, one can indeed completely recover their positions. 
In Figure 1.6, we show all tensors for the last example in more detail. Obviously, 
there are not only anisotropies; rather different morphological directions are rele- 
vant. The volume concentration around the center of mass is relatively high, since 
the volume is dominated by the two largest subclusters, being close to each other. 
It is quite the opposite with W 2 ' '■ P2 is relatively far away from all the clumps, 
the curvature around P2 is dispersed leading to a high value of tensor's trace. 
The main relevant direction connects the pair of large clumps with the third 
subcluster. The tensor W^' 1 is always isotropic, as shown in Equation ( 1.13 ), its 
trace is connected to the volume. Obviously, W 2 ' and Wf' 2 complement each 
other; indeed they add up to a diagonal tensor (recall the linear relationships 
described in |jt3|). Their eigenvectors' alignment with the coordinate axes is due 
to the fact that all the subclumps are aligned parallel to the axes as well. Oth- 
erwise, W 2 L,L and Wi would average over the orientations of the subclumps. 
Finally, W® ' 2 is isotropic, its trace equals the Euler characteristic. 
The detailed description of patterns given by the Minkowski valuations can be 
used in practical applications, as we shall see in the next section. 




Fig. 1.6. We show all of the tensors W 2 ' , W 2 ' , Wl' 1 , and W°' 2 for v = 1,2 - 
they span a base of the second-rank tensor space. The tensors Wl' a were divided by 
Wv'° (i — 0, ..,2); the tenors with a twofold position weighting (W 2 '°) are calculated 
around the corresponding centroid. All of the other tensors are translation-invariant. 
Therefore we depict them outside the pattern on the right-hand-side. Their eigenvalues 
are multiplied by a factor of 6 in order to yield a more convenient representation. Wq'°: 
solid line; W?>° ', W^ 1 , W°' 2 : dashed line; W 2 2 ' , W 2 M , W°' 2 : dotted line. 
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1.5 Physical applications of higher-rank Minkowski 
valuations 



1.5.1 The inner structure of spiral galaxies 



The classification of spiral or disc galaxies into specific types (the Hubble types, 
e.g., |]l9| ) is an important step towards their physical understanding; there are 
competing theoretical pictures trying to explain the formation of spiral arms 
[ p3|JT5[ |. In Figure 1.7 we display a real galaxy observed face-on together with a 
simulation according to a percolation model |H}] . So far one has mainly relied 
on an assignment by eye in order to classify observed disc galaxies and to com- 
pare observational data with theoretical models |2£| . One of the reasons for this 
is that spiral arms are easily detected by the human eye, but are rather hard 
to quantify. Here we use Minkowski valuations to reveal aspects of the inner 



structure of spiral galaxies (see 36, il] , p2[ for some discussion of how to quantify 
morphometric features of spiral galaxies, and especially |^,^|). 
In order to render an optical galaxy image accessible to a Minkowski analy- 
sis, it is helpful to understand it as a two-dimensional surface brightness field 
u(x). We may then estimate the Minkowski valuations of the excursion sets 
M v = {x|it(x) > v}. Contrary to large scale cosmological data (e.g. the distri- 
bution of galaxies or clusters of galaxies), such a density field is neither homo- 
geneous nor isotropic; rather it reveals a point of considerable symmetry: the 
center of mass. We expect the curvature centroids to coincide within this point 
and therefore concentrate on second-rank Minkowski tensors. 
In order to detect a possible overall anisotropy, we divide the smaller eigenvalue 
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tion, e = 1 — X can serve as a measure of the anisotropy. Such an anisotropy 
may reflect a physical anisotropy, induced, for instance, by the presence of two 
massive spiral arms, but may also indicate that the galaxy is not observed face 
on, but rather under an inclination angle. 

As the tensors W$'° (v = 0,..,2) weight the volume or the surface, respectively, 
with xx, their eigenvalues are increased, if the configuration is more widespread 



(see Fig. 1.8). This fact can be used to quantify the concentration of a density 
field around its center of mass, by calculating the trace of W^'°. Therefore an- 




Fig. 1.8. The tensors Wq'° (solid line) and W±'° (dashed line) for a wide spread (right 
panel) and a compact (left panel) configuration of circles. The eigenvalues depend 
quadratically on the separation of the circles. 



other characteristic feature of disc galaxies, namely the way in which the radial 
surface brightness distribution declines (with increasing radius), can be analyzed 
by scanning through various threshold values v. We applied (suitably normal- 
ized) concentration parameters P v (built from Tv(W^' )) to the images of both 
real galaxies (C. Moellenhoff, private communication) as well as to galaxies pro- 
duced by a percolation model (0); see jl3) for details). These morphological 
descriptors reveal strong luminosity concentration for real gala xies, which could 
not be reproduced by any of the model galaxies (see Fig.O) with reasonable 
choices for the model parameters. Using the morphology of observed disc galax- 
ies, one can therefore rule out the percolation model (at least in its simple form) 
as a viable attempt to explain the formation of spiral arms. 



1.5.2 The morphology of galaxy clusters 



As mentioned in Section 1.2, the morphology of clusters may constrain the values 
of the cosmological parameters. Intuitively, one may expect clusters in cosmologi- 
cal low-density models to have less substructure than in high-density models |3^| . 
Using Minkowski valuations, we check whether this effect can be observed for 
the cosmological GIF simulations M. We analyze the clusters' projected mass 



distributions such as shown in Figure 1.2. Observations from the gravitational 
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Fig. 1.9. Results of a quantitative morphometry for spiral galaxies using Minko wsk i 
tensors. We consider the galaxy NGC 1232 (shown in the left panel of Figure 1.7) 
and a sample of model galaxies following a percolati on m odel (an example of a model 
prediction is displayed in the right panel of Figure 1.7). We plot two concentration 
parameters (as described in the main text) versus decreasing surface brightness thresh- 
old values, i.e. increasing surface area. The error bars indicate la fluctuations around 
the expectation value of the model, whereas the single line refers to NGC 1232. The 
parameter Ca reveals significant differences in the concentration between the surface 
brightness profiles of real and model galaxies. The particular percolation model used 
thus is incompatible with real data. 



lensing through clusters |3^JTj| trace approximately the same cluster informa- 
tion. 

In order to make the simulation data, which consist of "particles tracing mat- 
ter", accessible to the Minkowski valuations, we construct a pixelized density 
field u(x) and smooth it with a Gaussian kernel of widthQ A = 0.3/i _1 Mpc. The 
resulting field «a(x) contains the overall shape of the cluster. Again we estimate 
the Minkowski valuations of the excursion sets, which contain all pixels above 
a given density threshold v. Our method to estimate the Minkowski functionals 
from the grid is based upon Crofton's formulae pt||37P ] . For each cluster we av- 
erage a function of the Minkowski valuations / over different density thresholds: 

-— (/) = / du x f (1.16) 

in order to obtain a compact set of parameters. Here (u\) 1 and (u\) 2 denote 
appropriately chosen integration limits. We consider five types of substructure. 



4 Length scales are given in units of h~ 1 Mpc. One Megaparsec (Mpc) equals about 
3.26 million light years, the constant h ~ 0.65 accounts of the uncertainty in the 
measurements of the Hubble constant. 
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1. The dumpiness C = y ((—Wa — l) ) counts the number of subclumps in 
quantifying how much the Euler characteristic Wi / tt deviates from one. 

2. The construction of the asymmetry parameter A is based on the fact that 
curvature centroids which do not coincide with each other, indicate asym- 
metry. Therefore, we calculate the perimeter of the triangle constituted by 
the curvature centroids, 2Wi (A (pj)), A= (Wi (A (p»))). 

3. The ratio of the tensors' eigenvalues, 71 (W£' a ) /j2 (W„ ,a ) tells us to which 
extent anisotropics are present. Here we consider the tensor W-f '° and define 

( 71 (W?'°) /Ta(w?' )). 

4. In general the curvature centroids p„ wander in space if the density thresh- 
old is varied; their fluctuations therefore measure the shift of morphological 
properties. Focusing on pi we define S= ((pi — (pi)) 2 )- This parameter is 
related to the centroid variation used in cluster investigations before p7|. 

5. Likewise, the eigendirections (parametrized by an angle a) of the tensors are 
changing; this generates the twist of morphological properties. Restricting 
ourselves to W^'° we employ the parameter T= ((a — (a)) 2 ). 

For more technical details, see Q. For each time period (or cosmological red- 
shift) we estimate the order parameters for three projections per cluster and 
average over all clusters in one cosmological model in order to extract a typical 



morphology. In Figure 1.10, we display the evolution of the averaged morpho- 
logical parameters vs. redshift. One recognizes immediately a relaxation process, 
leading to less substructure (and lower values of the morphological order param- 
eters) as time goes by (from right to left). Furthermore, clear differences between 
the cosmological background models arise. Indeed the theoretical predictions ac- 
cording to which the low-density clusters (OCDM) exhibit less substructure, are 
confirmed for the past; for redshift z = 0, the differences disappear. 
All of our order parameters show a similar behavior. The reason probably is that 
we average over different types of relaxation processes. An investigation of single 
clusters should reveal differences between, e.g., a bimodal merger dominated by 
two colliding subclumps and a more isotropic mass infall. 



1.5.3 The geometry of the electric charge distribution in molecules 

An important implication of the atom hypothesis is that the structures of mole- 
cules are to be explained in terms of their constituent parts, the atoms. Heuristic 
principles based upon this assumption work quite well in chemistry. However, 
if one tries to justify this hypothesis in the framework of quantum mechanics, 
the notion of an individual atom seems to disappear, since the negative charge 
distributions of the atoms merge, blurring the electron clouds which define the 
single atoms. 

Bader [|| proposed an elegant method to recover the notion of atoms in molecules. 
Since, from the quantum mechanical point of view, the charge density profile g(r) 
is the only relevant quantity, its geometry and topology must contain the notion 




Fig. 1.10. The average morphological evolution of galaxy clusters from the GIF sim- 
ulations. Our morphological order parameters dumpiness C, asymmetry A, anisotropy 
X, shift of morphology S and twist of morphology T are shown versus time. The time is 
given in terms of the observable cosmological redshift. Redshift z — means the present 
day, higher redshifts lead back into past. In order to get typical results for two cos- 
mological models (low-density model with an open geometry, called open Cold Dark 
Matter, OCDM: dashed line; high-density model with a massive r neutrino, named 
rCDM: solid line) we average over about 15 clusters and 3 orthogonal projections per 
cluster within each model. As the results show, our morphological order parameters 
quantitatively feature the differences between the models which are in accordance with 
the theoretical predictions. For details on the cosmological Cold Dark Matter models 
see H 



of atoms. Using the critical points of the density distribution g, where Vp = 0, 
one can define atoms as local density maxima together with a surrounding do- 
main enclosed by zero-flux surfaces. This geometrical approach can be extended 
to chemical bonds and has been justified on physical grounds. Moreover, it pro- 
vides us with means to visualize molecules in a proper way using the charge 
density profile, see Figure |l.ll| for an illustration. A distinguished role is played 
by the surface V 2 g = which separates the zone where chemical reactions are 
likely to happen from the rest of space. 

Obviously, this is an area of research where physics and geometry strongly inter- 
act, and the applications of higher-rank Minkowski valuations may yield useful 
information on the structure of the chemical bonding. First, the morphology of 
the whole molecules can be simply described in terms of the Minkowski val- 
uations of the density field. The atoms, which Bader defines H, may not be 
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Fig. 1.11. The water molecule H2O. The bright lines denote the bonds connect- 
ing the H-atom with both of the O-atoms. The crosses mark the points where the 
zero-flux surfaces separating the atoms intersect with the bonds. Moreover, some iso- 
surfaces of the Laplacian of the charge density, Ao , are shown. The images are from 
fittp:/ /www.nas.nasa.gov/^creon/papcrs/mgms96/ (courtesy C. Levit). The panels 
show different projections and focus on different iso-surfaces. 



compact in every case, so the single atoms can not be understood in terms of 
the Minkowski valuations. However, the fact that the atoms' positions are spec- 
ified in Bader's theory gives us the chance of probing the molecular structure 
with the Boolean grain method where the atoms' positions are decorated using 
spheres (either monodisperse or with radii mirroring physical quantities). Note, 
that at this point one really needs tensors, since in this case the locations and 
the orientations of subsystems are of crucial importance; and this local informa- 
tion lies beyond the scope of the scalar measures. Bader's theory constructs also 
a graph of bonds, that may be analyzed with Minkowski valuations straight- 
fowardly. 

But given the fact that both the molecules and the constituting atoms are de- 
fined in terms of geometry, one may go one step further and try to discover 
phenomenological principles allowing one to predict the shape of molecules given 
their constituents. 

These endeavors clearly are beyond the scope of this article. However, in order 
to illustrate to some extent how useful the Minkowski valuations may prove in 
this area, we consider briefly some of the Minkowski tensors for the H^ ion. 
In order to treat the H^ ion consisting of two protons and an electron in the 
framework of quantum mechanics, one often considers the protons as fixed in 
space (Born-Oppenheimer approximation) thereby reducing the problem to an 
effective one-particle problem of an electron moving in the potential of two pos- 
itively charged nuclei. A reasonable ansatz for the quantal wave function of the 
electron - at least suitable to illustrate what is basically happening - is a super- 
position of quantum states where the electron is bound to either of the protons. 
To be more precise, let us consider the protons as being separated by their 
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Fig. 1.12. The morphology of the bonding and the anti-bonding states of the electron 
within the ion. We explore the charge distribution of the electron via the Minkowski 
tensors of excursion sets from the density field. In the left panel, we show the trace 
of W-j 2 ' as a function of the charge fraction q enclosed within the density contour. In 
the right panel, we consider the anisotropy of the charge distribution. It is quantified 
via An = (jW?'°) xx - {Wi'°) yy ) /Tr (W?'°). The anti-bonding state, shown as the 
dashed line is much more anisotropic than the bonding state which wraps both atoms 
with a smooth electron hull. 



equilibrium distance of approximately two Bohr radiQ Let us label the protons 
with 1 and 2 and let ipi denote the unperturbed (well-known) hydrogen ground 
state wave function of an electron bound to proton i for i — 1,2. For symmetry 
reasons, only the linear combinations 

ip a = (ipi - ip 2 ) and ip b = —= (tp! + ip 2 ) (1-17) 

are of interest. Since ipb is favored on energetic grounds, it is a bonding state; 
consequently ip a is called antibonding state. 

The morphology of the electron charge distribution is expected to distinguish 
between both states. Let us consider excursion sets of the three-dimensional elec- 
tron charge density profile g = \ip\ 2 for both states. In Figure 1.121 we display the 

r2,0 ~ -' 



surface tensor W x ' as a function of the charge fraction enclosed within the ex- 
cursion set. We start from the very peaks of the charge distribution (high charge 
density thresholds) and move on to lower charge density thresholds enclosing a 
significant part of the whole electron charge. The trace of the tensor W^'° shows 
that the iso-surfaces of the bonding state are slightly more extended than those 

5 In the following we consider spatial distances in units of Bohr radii. The Bohr radius 
ao marks the typical scale of the hydrogen ground state wave function. For an intro- 
duction into the physical treatment of the ion one may consult, e.g. Complement 



Gxi in jTJj. 
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of the anti-bonding state. But on the whole, the differences are small: the con- 
centration of the charge distribution is basically determined by the distance of 
the protons. If we ask for morphologically relevant directions, the differences be- 
tween both states are much more perspicuous. The anisotropy of the azimuthal 
symmetric quantum states can be quantified via 

An=((wf°) ^/Tr(w?°) , (1.18) 

the protons lying at (x,y, z) = (—1,0,0) and (1,0,0), respectively. For small 
charge fractions, the tensor VKj 2 ' is dominated by the peaks around the pos- 
itively charged protons yielding a maximum anisotropy, but no essential dif- 
ference between both states. However, as we move to lower density thresholds 
and higher charge fractions, the bonding state isotropises much more quickly 
than the anti-bonding one. In the bonding state, therefore, the electron forms 
a cloud spreading out over both atoms in a quite isotropic way, whereas in the 
anti-bonding state, the signature of the two distinct protons remains still visible. 
The chemical bond therefore is mirrored by the geometry of the electron density 
distribution, as quantified via Minkowski tensors. 

This example is but an elementary first illustration of the nature of the chemical 
bond. In the spirit of Badcr's theory, one can possibly go much further in order 
to consider single atoms in molecules. Additivity can serve as a guiding princi- 
ple, since Bader found that there are stable atom configurations which persist 
through different types of molecules, presenting a sort of rigid module movable 
as a whole; this fact is well-known from chemistry, of course. 

1.6 Density functional theory based on Minkowski 
valuations 

Having shown the morphometric versatility of MVs, we conclude our discussion 
of applications by considering systems where integral geometric quantities en- 
ter a physical approximation. The quantities that come into play at this point 
generalize the Minkowski valuations further and may indicate a new direction of 
interest. 

Density functional theory became the standard method to study inhomogeneous 
classical fluids in the last 20 years. The main feature of such physical systems is 
the spatial variation of the equilibrium density p e q{~x) which minimizes the free 
energy fl. Density functional theory is based on the exact result that the free 
energy of the inhomogeneous fluid can be expressed as a functional i?[p(x)] of 
the density p(x). Thus, all relevant thermodynamic properties such as surface 
tensions or phase transitions of confined fluids can be calculated from i7[p e9 (x)]; 
even equilibrium correlation functions that describe the microscopic structure of 
inhomogeneous fluids can be determined by derivatives of the functional. 
In general, the form of the functional J?[p(x)] is not known; deriving the exact 
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functional would be equivalent to solving the statistical mechanics for fluid sys- 
tems exactly which seems to be an impossible task. Therefore, one has to restrict 
oneself to reasonable approximations for i7[p(x)]. For convenience, one separates 
an excess intrinsic free energy functional f ea; [{/5i}] from exactly known external 
contributions 

fl[p(x)] = y~dx (f ex [p(x)] + fc B Tp(x) (lnyl 3 p(x) - 1) + p(x) (K(x) - /;)) 

(1.19) 

where V(x) is an external potential, /i the chemical potential, (3 = l/(fesT) the 
inverse temperature T, and A the thermal wavelength of the fluid, ks denotes 
the Boltzmann constant. 

The most elaborated and reliable functional for hard-sphere fluids is that devel- 



oped by Rosenfeld 1 33 based on the local scalar and vector- valued Minkowski 
functionals 

W»'«[p(x)] = -L f dS^s^t (k x , .., Kd _r) x*>n* p(x) . (1.20) 



In contrast to the definition (1.6) and the applications considered in the previous 
sections, the functionals here are locally weighted with the density p(x). Notice, 
that for inhomogeneous densities the vector valuations W^ ,:L [p(x)] in general do 
not vanish for a convex body K and that the functionals do inherently depend 
on the location of K . The following ansatz for the excess free energy functional 
proves to be useful: 

/3Up(x)] = A«<°) Wf° + / 2 «°) <%r + /3«°) «°) 3 
+/ 4 «'°) wf- V, - 1 + / 5 «'°) <-°«'V 2 ' 1 ) , 

(1.21) 

where we drop the argument [p] for convenience. This ansatz is inspired by the 
observation that the physical dimensions are [W 3 ,0 ] = [W x ' ? W 2 ' 9 ] — [{W® ' 9 ) 3 ] = 
\Pfex] — [length] ~ 3 , that [W® ,0 ] = 1 is a dimensionless weighted density. 
To obtain the functions fi(W : °) one may apply an argument from scaled particle 



theory for a homogeneous fluid |pl| , namely, that the excess chemical potential 
for big spheres of radius a/2 — ► oo is Yaa.fi /V — > p with the volume V = §c 3 , 
the chemical potential fi = and the pressure p. Within the present ansatz 



(1.21) one finds in this limit the reversible work to create a cavity of volume 

df ea 
aw" 



V, fi — rZ^u.u q , and therefore in the limit a/2 — > oo the chemical 



potential ji/V —> ?.l'5,<> ■ Using the equation of state fl — ~pV, i.e., p = — f, 



ex 



fw" ^3 ' one S e t s t ne wen known 'scaled particle' differential equation 
dw o,o ^ dW^° dW?' 1 1 dW° A 2 3 

(1.22) 



1 Vector- and tensor-valued descriptors for spatial patterns 



21 



Within the ansatz (1.21) the solution of the differential equation reads 



/3/e*[p(x)] = - 



3W. 



o.o 



4-7T 
7-0,0 



ln(l -<•") + 



9 W?'°W°'° 



47T 



i - w; 







0,0x2 



+ 



9W^ (W 1 U,U ) 



(1.23) 



S7T 



(i-wy- u ) 



0,0n 2 



This functional yields excellent results for thermodynamic quantities and for 
structu res i n the fluid phase of hard spherical molecules j}5|. For instance, in 
Figure 1.13 the inhomogeneous density profile p{z) close to a hard wall is shown 
as compared to simulation results. 



Unfortunately, the functional in Eq. (1.23) was derived only for hard spheres in 
three dimensions; however, a proposal has been made for general convex bodies 




[ p4| based on an approximative decomposition of f ex in the low density limit. 
Here, we show how Rosenfeld's idea of decomposing f ex for spheres j33[ may be 
extended exactly to non-spherical convex bodies by applying integral geometric 
techniques such as kinematic formulae and an explicit local expression for the 
Eulcr characteristic \{K n K') of two overlapping bodies pq] . The main idea 
is to decompose the exact expression for the local excess free energy f ex [p(x.)] 
at low densities in terms of the vector-valued Minkowski functionals. To this 
end one has to evaluate the configurational integral J dK = J dR J dx over all 
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orientations R and positions x of a body K = RK X , namely, 

y J dx/?/ ex [p(x)] -> y / 1 dX' X (A' x n^ x p(x)p(x') 

= / dx <°[p(x)]<°[p(x)] 

dx 3<'°[p(x)]Ty 2 '°[p(x)] - 3<> 1 [p(x)]W 2 °' 1 [p(x)] 

, /tx - k 2 (vivi - v 2 v 2 ) i ■ (n 'n % 



dx / dX x p(r) / d<p(r')' 



1 + nn' 

(1.24) 



where Vj denotes the principal direction of the curvature Ki at r and n (n') 
are the normal vectors of RK X (R'K X ) at r (r'), respectively. The integral over 
rotations R and over the surface of the body K x centered at x is written as 
J r d-f^x = / di? J re g RK dS''' -1 . For spheres with k,\ — k 2 one recovers immedi- 



ately the low-density limit of the intrinsic free energy given by Eq. (1.23). But 



unfortunately the last term on the right hand side of Eq. (1.24) vanishes only for 
spheres, so that only for hard spheres a finite number of Minkowski functionals 
is, by accident, sufficient to decompose the excess free energy into contributions 
stemming alone from the individual bodies, whereas for non-spherical shapes in 
arbitrary dimensions an infinite number of fundamental curvature measures is 



required. Nevertheless, based on the most general analytical expression (1.24), 
several approximations by a finite number of tensor- valued Minkowski function- 
als (including the proposal in [|34"|) may be proposed and numerically tested 
which make Rosenfeld's density functional applicable to general convex bodies, 
in particular, to spherocylinders and ellipsoids p6[. 



1.7 Conclusions 



Concluding, we argue that a new perspective emerges as regards the Minkowski 
functionals. Instead of marking a small number of distinguished measures (which, 
of course, they do in a certain sense), the Minkowski functionals are embedded 
into a more extensive class of Minkowski valuations. These measures generalize 
the Minkowski functionals without leaving the framework of integral geometry. 
The essential step to extend the scalar functionals is to replace the motion invari- 
ancc by motion covariance. This step yields a hierarchy of tensor- valued measures 
with the Minkowski functionals at the bottom. The higher-rank valuations are 
moments of the Minkowski functionals where local geometry is weighted with 
either the position or the normal vector. The full, comprehensive description of 
the tensor-valued measures is still a mathematical challenge, although some im- 
portant results concerning the integral geometry of tensor measures have already 
been achieved, especially in low dimensions (see p4|). 
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The extensions of the Minkowski functionals are not only mathematically inter- 
esting; rather they can be motivated on physical grounds. Since the Minkowski 
functionals proved to be useful for a number of applications, we are justified 
to expect that also their generalizations may be suitable tools. Using physical 
applications at various length scales we exemplified how higher-rank Minkowski 
valuations complement the information featured by the scalar functionals. From 
this point of view we hope that the higher-rank Minkowski valuations will become 
a standard tool for morphometry and some other fields of research concerned 
with pattern formation. 
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